Deconfinement in Matrix Models about the Gross— Witten Point 
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We study the deconfining phase transition in SU(N) gauge theories at nonzero temperature using 
a matrix model of Polyakov loops. The most general effective action, including all terms up to 
two spatial derivatives, is presented. At large N, the action is dominated by the loop potential: 
following Aharony et al, we show how the Gross- Witten model represents an ultra-critical point in 
this potential. Although masses vanish at the Gross- Witten point, the transition is of first order, as 
the fundamental loop jumps only halfway to its perturbative value. Comparing numerical analysis 
of the A = 3 matrix model to lattice simulations, for three colors the deconfining transition appears 
to be near the Gross- Witten point. To see if this persists for N > 4, we suggest measuring within 
a window ~ 1/A 2 of the transition temperature. 

PACS numbers: 



I. INTRODUCTION 

As an SU(N) gauge theory is heated, it undergoes a 
phase transition from a confined, to a deconfined, phase. 
The standard order parameter for deconfinement is the 
Polyakov loop in the fundamental representation |]J . On 
the lattice, numerical simulations universally measure the 
expectation value of this loop to see where the deconfin- 
ing transition occurs |2|. This is only of use at finite 
lattice spacing, though, since in the continuum limit, it 
appears that the expectation values of all bare loops van- 
ish 0. 

Recently, two techniques have been developed to mea- 
sure renormalized loops, which are nonzero in the con- 
tinuum limit LjJI - A Polyakov loop represents the prop- 
agation of an infinitely massive, "test" quark. On a lat- 
tice, even an infinitely massive field undergoes an additive 
mass shift, which diverges in the continuum limit. This 
mass shift generates a renormalization constant for the 
loop: as for any bare quantity, this must be divided out in 
order to obtain the expectation value of the renormalized 
loop. 

Following Philipsen 5] , one method computes the sin- 
glet potential 0, which is equivalent to computing the 
Wilson loop at nonzero temperature. When the Wilson 
loop is narrow in the spatial direction, it can be computed 
in perturbation theory, allowing for the renormalization 
constant of the loop to be extracted [(J Q • 

The second method involves the direct measurement 
of one point functions of Polyakov loops. The drawback 
is that these must be measured on several lattices: all 



at the same physical temperature, but with different val- 
ues of the lattice spacing Q. While these two methods 
are rather distinct, for the triplet loop in a SU(3) gauge 
theory, they agree to within w 10% over temperatures 
from Td — > 3Td, where Td is the transition temperature 
for deconfinement |(J. 

In an asymptotically free theory, the expectation value 
of a Polaykov loop approaches one at high temperature. 
It is perfectly conceivable that when a gauge theory de- 
confines, that it goes directly to a gluon plasma which 
is close to perturbative, even at temperatures just above 
T d . (An example of this is Fig. © of Sec. ITITdT ) In- 
deed, if the deconfining phase transition is strongly first 
order, this is what one might naively expect. From lattice 
simulations, while the deconfining transition is of second 
order for two colors 

urn, 

it is of first order for three col- 
ors 0, El , and becomes more strongly first order as the 
number of colors increases. Results for four 0, [T2I Il3| , 
six, and eight colors [T3 . Il3| suggest that at large N, the 
latent heat is proportional to an obvious factor of ~ A^ 2 , 
from the overall number of gluons in a deconfined phase, 
times a constant. 

For three colors, the lattice measurements of the renor- 
malized triplet loop are admittedly preliminary. With 
that caveat, the expectation value of the renormalized 
triplet loop is found to be near one, ~ .9 ± 10%, for tem- 
peratures which are as low at ~ 3Td- This agrees with 
resummations of perturbation theory, which work down 
from infinite temperature, but consistently fail at a tem- 
perature which is something like ss 3 — ► 5Tj [l4|. 

Just above the transition, however, the renormalized 
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triplet loop is far from one: it is only w 0.4 ± 10% at T^~. 
This suggests the following picture. For temperatures 
above « 3Xd, the deconfined plasma is, after resumma- 
tion, nearly perturbative. For temperatures between 
and w 3Xd, however, the theory is in a phase which is 
deconfined but non-perturbative. The existence of such 
a non-perturbative (quark) gluon plasma may be indi- 
cated by the experimental data on heavy ion collisions at 
RHIC jll. 

At large N, the deconfining transition has been studied 
using string and duality methods [161 Il7l Il8l lltl |2(I |21| . 
At infinite N, by restricting the theory to live on a space- 
time where space is a very small sphere, it is possi- 
ble to compute the Hagedorn temperature analytically 
Ellil El 111 EHHl. In this limit, Aharony et al. [3 
catagorized the possible relationships between the decon- 
fining, and the Hagedorn, temperatures. 

In this paper we develop a general approach to the 
deconfining transition, based upon a matrix model of 
Polyakov loops. We begin in Sec. [H] with the most gen- 
eral effective Lagrangian, including all terms up to two 
spatial derivatives. At first sight, matrix models of loops 
appear to be an ine lega nt t ype of nonlinea r sig ma model 

[2a HE E3, 12E m, [23, [23, [23, 113, HH H3, Uli . This 

is not a matter of choice, but is dictated by the physics. 
For ordinary nonlinear sigma models, in the chiral limit 
there are no potential terms, only terms involving spa- 
tial derivatives. These can be studied in a mean field 
approximation, but this is only valid when the number 
of spacetime dimensions is large [^[^El^. ^ n con trast, 
the effective action for loops starts with a potential, inde- 
pendent of spatial derivatives. Because of this, a system- 
atic large N expansion can be developed in any number 
of spacetime dimensions. 

The large N expansion of the loop potential is devel- 
oped in Sec. IIIII We show how the Vandcrmonde de- 
terminant, which appears in the measure of the matrix 
model, gives a type of potential term 0, 13H l39| . also 
contributes a potential term. At infinite N, the vacua 
of the theory are then given by the stationary points of 
an effective potential, which is the sum of the loop, and 
what we call the Vandermonde, potentials. 

In the loop potential, the adjoint loop is a mass term, 
while loops in higher representations represent interac- 
tion terms. We show how the Gross-Witten model 
[3^. Eof. which arose in a very different context, corre- 
sponds to a loop potential which is purely a mass term, 
with no higher terms. We term the point at which the de- 
confining transition occurs as the "Gross-Witten point" 
0. At this point, the expectation value for the loop 
in the fundamental re pres entation, £, jumps from to 
precisely \ 0, H S H3- The loop potential is most 
unusual, however: it vanishes identically for i between 
and 5, and is only nonzero when t > | (see Fig. ^ be- 
low). This is only possible because of the contribution of 
the Vandermonde potential. This potential, which is not 
analytic in I at I = i, produces a deconfining transition 
which is thermodynamically of first order d EH!!. This 



flat potential also implies that at the transition point, 
masses vanish, asymmetrically, in both phases (these are 
the masses for the fundamental loop, along the direction 
of the condensate) 3]. If a background field is added, at 
the Gross-Witten point the first order transition turns 
into one of third order, no matter how small the back- 
ground field is mm. 

Aharony et al. analyzed the transition when quartic 
interactions are included 01 ■ They suggested that in the 
space of loop potentials, the Gross-Witten model repre- 
sents a tri-critical point. We refine their analysis, and 
extend it to include arbitrary interactions of the funda- 
mental loop. We discuss why the Gross-Witten point 
is more accurately described as an ultra-critical point, 
where all higher interactions vanish. Physically, this is 
defined as the unique point where the transition is of first 
order, and yet masses vanish; it only arises at infinite N. 

We analyze the matrix integral by introducing a con- 
straint field for the fundamental loop |35j. Various con- 
straint fields can then be integrated out in different order. 
One order is equivalent to computing the Vandermonde 
potential by Legendre transformation |18l Il9l |2(1 l21j . 
We show that the other order is equivalent to what was 
called, previously, a mean field approximation pj. l4ll E3 , 
m. If the potential includes just the adjoint loop, we 
check that while these two methods give different effec- 
tive potentials, that at any stationary point, both the 
expectation values of the loop, and masses, agree. 

For three colors, the surprising aspect of the lattice 
data is that the value of the renormalized triplet loop 
at T£, w 0.4, is close to |. Further, masses associated 
with the triplet loop — especially the string tension — 
decrease significantly near the transition |10| . This sug- 
gests that the transition for three colors is close to the 
Gross-Witten point at infinite N. 

The lattice results provide further evidence for the util- 
ity of using a large N expansion at N = 3. The expecta- 
tion values of renormalized sextet and octet loops, which 
are expected to be ~ 1/N or smaller, never exceed w 25% 
at any temperature, and are usually much less 

In Sec. IIVI we study the loop potential for finite TV. 
Following Damgaard et al. 0] , we start by showing that 
when fluctuations in the matrix model are neglected, 
that the expectation value of any loop — including those 
which arc Z(N) neutral — vanish in the confined phase. 
This is a striking difference between matrix models and 
more general Z{N) symmetric theories pHj . 

The matrix model is merely a two dimensional integral, 
which can be studied numerically. In Sec. IIVI we inves- 
tigate theories close to the N = 3 analogy of the Gross- 
Witten point, including cubic interactions. The matrix 
model gives the expectation values of loops in arbitrary 
representations. We also compute the mass of the loop, 
which is related to the (gauge invariant) Debye mass, as 
obtained from the two point function of Polyakov loops. 
We discuss how measurements of both the Debye mass, 
and the renormalized triplet loop, help probe the loop 
potential. 
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We also use the matrix model to compute difference 
loops for the sextet, octet, and decuplet loops. If only 
the loop potential is included, the difference loops from 
lattice simulations are not well reproduced. Presum- 
ably it is necessary to include fluctuations, due to kinetic 
terms, in the effective theory. Nevertheless, up to these 
small corrections, < 20% at all temperatures, it is clear 
that the deconfining transition for three colors is close to 
the Gross- Witten point at infinite N. 

In the Conclusions, Sec. [3 we discuss why lattice data 
[Til IT2I Il3| may indicate that the deconfining transitions 
for N > 4 are also near the Gross- Witten point. We 
suggest that one must be very close to the transition to 
see this, when the reduced temperature \T — Td\/Td ~ 
1/N 2 . 



II. EFFECTIVE MODELS 

A. Nonlinear Sigma Models 

We begin by briefly reviewing effective theories for non- 
linear sigma models j^. For defmiteness, consider the 
model appropriate to chiral symmetry breaking for Nf 
flavors, with a global symmetry group of |46( 



G f = SU L (N f ) x SU R (N f ) 
and a field U, transforming as 



(1) 



(2) 



where £Il,r are SU L,Fi(Nf) transformations. Because the 
left and right handed chiral rotations are distinct, in the 
chiral limit, terms such as the trace of U, tr U, cannot 
arise. Thus in the chiral limit, there is no potential for 
U, and there are only derivative terms. These start at 
second order: 

C = fl tr \d»U\ 2 + . . . ; (3) 

subject to the constraint that U is a unitary matrix, 

U^U = l Nf . (4) 

The series then continues with terms of quartic order in 
derivatives pfij ]. 

After chiral symmetry breaking, what remains is a vec- 
tor symmetry of SUv{Nf), with Ql = = ^v- Then 
there is a potential for U possible, 



V = mi tr U 



(5) 



This potential is manifestly proportional to the pion mass 
squared, since it must vanish in the chiral limit. Near the 
chiral limit, the kinetic terms in @ dominate over the 
potential term in (|5j) if the volume of space (or space- 
time) is large. If the volume is small, however, the po- 
tential term can dominate , as the modes with nonzero 
momentum are frozen out |47| . 



Other sigma models are constructed by allowing the 
transformations and to be equal. For example, 
take U to be a SU(2Nf) matrix, and impose a further 
constraint, such as that the trace of U vanishes, tr U = 0. 
In this case, the symmetry is that of a symmetric space, 
Gf = SU(2Nf)/S(U(N f )xU(N f )) [MEI. Because the 
trace of U is constrained to be a fixed number, there is 
no potential possible, and the action is given entirely by 
kinetic terms, as in J3J). If the trace is constrained to be 
some other value, then the symmetry group changes, but 
there is still no potential possible. 

These nonlinear sigma models are renormalizable in 
an expansion about two space-time dimensions. There 
is a phase with broken symmetry above two dimensions, 
where the expectation value of U is nonzero. This ex- 
pectation value is generally not proportional to the unit 
matrix, and there are Goldstone bosons in the broken 
phase. In the chiral limit, as there is no potential for U, 
it is not easy studying the possible patterns of symmetry 
breaking directly in the nonlinear form of the model (lin- 
ear models are usually more useful) . Symmetry breaking 
can be studied using mean field theory [35|,|3(|. To do so, 
the continuum form of the theory is replaced by a lattice 
form, with matrices U on each site i. In that case, the 
next to nearest neighbor interaction at a site i, for unit 
lattice vector n, becomes, in mean field approximation, 



tr {UiUi+n) -> tr {Ui (U)) 



(6) 



In this case, mean field theory is only applicable when the 
number of nearest neighbors, or space-time dimensions, 
is very large [3{| |3(j . 



B. Polyakov Loops: Preliminaries 

We discuss the deconfining transition at a nonzero tem- 
perature T, for a space of infinite volume. The analysis 
is more general, though, and can be easily extended to 
the case in which space is of finite extent, etc. 

The thermal Wilson line is 



Ltc(x) = V exp ig 



l/T 



(7) 



We follow our previous notations and conventions . We 
define the thermal Wilson line at a point x in space, let- 
ting it run all of the way around in imaginary time, from 
to r = l/T. Otherwise, V denotes path ordering, g 
is the gauge coupling constant, and Aq the vector po- 
tential in the time direction. The are the generators 
of SU(N) in a representation 1Z, which is taken to be 
irreducible. 

The Wilson line transforms homogeneously under local 
gauge transformations, 



L n (x) -> rt n (x, l/T) L n (x) n n {x,0) 



(8) 



For gauge transformations which are periodic in imag- 
inary time, we form a gauge invariant quantity, the 
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Polyakov loop in 1Z, by taking the trace, 

&n = -7— tr Lk . (9) 

&n 

We define the Polyakov loop as the normalized trace of 
Lrji, with dfi equal to the dimension of 1Z. The advan- 
tage of this is that in an asymptotically free theory, all 
Polyakov loops are of unit magnitude in the limit of in- 
finitely high temperature. We denote representations by 
their dimensionality, with two exceptions, which differ 
from previous use [3j . The Wilson line in the fundamen- 
tal representation is L, with the fundamental loop 

L N = 1 tr L . (10) 

The adjoint loop is denoted 

^^-^(ItrLl 2 -!) . (11) 

Besides local gauge transformations, which are strictly 
periodic in time, in a pure gauge theory there are 
also global gauge transformations, which are only pe- 
riodic up to an element of the center of the gauge 
group. For an SU(N) gauge group, the center is Z(N) . 
In the fundamental representation, the simplest global 
Z(N) transformation is 

fi(f,l/T) = e 2m ' N n(S,0) . (12) 

Defining the charge of the fundamental representation to 
be one, the charge of an arbitrary representation, e-jz, 
follows from its transformation under 1)12(1. 

in - e 2 ^ /N in ■ (13) 

As a cyclic group, the charge e-n is only defined mod- 
ulo N. A special set of representations are those with 
zero Z{N) charge, which we denote as H®. The simplest 
example is the adjoint representation. 

Below the deconfining transition temperature T d , the 
expectation values of all loops with nonzero Z(N) charge 
vanish, 

(M=0 , T<T d , K^Ko, (14) 

and the theory is in a Z(N) symmetric phase. 

The global Z{N) symmetry is broken in the deconfined 
phase, as loops in all representations condense, 

{In) ^0 , T>T d , V U . (15) 

In the limit of large N, factorization 

[HUES! implies 

that all expectation values are powers of that for the 
fundamental, and anti-fundamental, loops: 

(e R } = {tN) p+ {t N y- , n = 00. (i6) 

The integers p+ andp_ are defined in 0. For the adjoint, 
for example, p+ = p- = 1; in general, the Z{N) charge 



e-jz = p + — P-. Factorization is exact only at infinite N, 
but holds for all temperatures. Corrections to factoriza- 
tion are of order ~ 1/N at large N; for the adjoint, they 
are~ 1/iV 2 0. 

In general, loops with zero Z(N) charge, such as the 
adjoint, can have nonzero expectation values at all tem- 
peratures, including in the confined phase. At infinite N, 
factorization implies that they vanish below T dl because 
the expectation value of the fundamental loop does. An 
important question to which we shall return is the extent 
to which Z(N) neutral loops condense in the confined 
phase. 



C. Matrix Models of Polyakov Loops 

In the same spirit as familiar for effective models of chi- 
ral symmetry breaking 0,E3> discussed in Sec. Ill Al we 
construct an effective theory for deconfincmcnt, applica- 
ble over distances greater than the inverse temperature. 
Thus fields depend only upon the spatial coordinates, x. 
We take as the basic field the Wilson line, L(af); as that 
is a gauge covariant field, we need also include the spatial 
components of the vector potential, Ai (x) . 

We are then led to construct the most general ac- 
tion consonant with the relevant symmetries, which are 
those of global Z(N) , and local SU(N), transforma- 
tions. We will allow the vacuum to spontaneously break 
the Z{N) symmetry, but do not allow SU(N) to break. 
(It would be interesting to extend the following analysis 
to include a Higgs phase, as is appropriate for symmetry 
restoration in the electroweak interactions.) 

We start with potential terms, independent of any 
derivatives in space-time. Because of the left-right sym- 
metry, for sigma models there is no potential in the chiral 
limit. In constrast, loop models always have a potential. 
To be invariant under SU(N) transformations, terms in 
the potential must be a sum of loops; to be invariant 
under Z(N) , these must be Z(N) neutral: 

V(L) - -m 2 l ad +J2 Ks~Re£ s ; (17) 
Sen 

Loops can be complex valued when N > 3, so in the 
potential we take the real part of each loop. 

By using group theory, this potential can be rewritten 
in many different ways. As higher powers of any loop can 
always be reexpressed as a linear sum over loops in other 
representations, this potential is the most general form 
possible. 

In V(L), the adjoint loop looks like a mass term for 
the fundamental loop, while loops in higher representa- 
tions look like interactions of the fundamental loop (plus 
new terms, such as tr L 2 /A, etc.). The loop potential 
differs from that of ordinary scalar fields, though. In 
the potential of a scalar field, the coupling constant for 
the highest power of the field must be positive in order 
for a ground state to exist. This isn't true for the loop 
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potential, as each loop is bounded by one: there is no 
constraint whatsoever on the signs of the mass squared, 
nor on any of the coupling constants. For example, we 
choose the sign of the adjoint loop to have a negative 
sign (so that m 2 ~ T — Td, as for ordinary spin systems). 
Then m 2 —>■ — oo drives one to the confined phase, while 
m 2 — > +00 drives one into the deconfined phase. Further, 
it is consistent to have just a mass term, with no other 
terms in the potential, setting all Kg = 0. Thus while the 
adjoint loop looks like a mass term for the fundamental 
loop, because the basic variable is L, and not i^, the 
theory is still nontrivial with no "interactions" , Ks = 0. 

At each point in space, the Wilson line can be diago- 
nalized by a local gauge transformation, and so depends 
upon N — 1 eigenvalues. Thus instead of L, it is also 
possible to chose a set of N — 1 loops, and rewrite the 
loop potential in terms of these. We find this of use in 
Sec. IV, (J70J). 

The loop potential can be computed analytically in 
two limits: when the temperature is very high |33l l34l . 
and when space is a small sphere 0, 0, lisl Il9i I20I 
I21I ]. For both cases there is a large mass scale, either 
the temperature, or the inverse radius of the sphere, so 
that the effective gauge coupling is small by asymptotic 
freedom. In perturbation theory, the loop potential is 
computed from the one loop determinant in the presence 
of a b ackgrou nd gauge potential A a ~ log L fl7L fla . 
IT9ll20ll2lLl3ll3i|. At very high temperature (and infinite 
spatial volume), while the result can be written in terms 
of an infinite polynomial of (traces of) powers of L, as in 
p7JI. the sum can be explicitly performed. The result is 
just a a simple quartic potential of Aq |33|. 

In a perturbative regime, it is natural that the loop 
potential is more transparent in terms of elements of the 
Lie algebra, the A , instead of elements of the Lie group, 
the L. This suggests that an effective theory of loops is 
not especially convenient when the theory is essentially 
perturbative. For SU(3), resummations of perturbation 
theory suggest that the perturbative regime appears for 
T w 3T<j |9- Thus for SU(3), loops are only useful for 
temperatures below w 3T c i- 

The loop po tential can also be computed on a small 
sphere [H M, 13 H3, E| . There are now two scales in 
the problem — the radius of the system, and the tem- 
perature — so the result is much more complicated than 
for infinite volume. The result is a sum like that of (|17|l . 

Kinetic terms, involving two derivatives, are con- 
structed similarly. While sigma models only have one 
term with two derivatives, due to the change in sym- 
metry, loop models have an abundance. One class of 
terms involves covariant derivatives of the Wilson line: 

£i(L) = 4 tr l^ L l 2 ( 1+ J2 X s Rci s) (18) 

5 V SeKo ) 

The first term is a sort of "electric loop" : it is the original 
electric field of the gauge theory, rewritten in terms of 
loops. At tree level, the coupling constant g equals the 



gauge coupling g, but in general, it is an independent 
coupling constant of the effective theory. Besides the 
electric loop, there is an infinite series of terms involving 
loops in Z(N) neutral representations. The couplings 
of magnetic fields is similar to 118|) : there is an infinite 
series of Z(N) neutral loops which couple to the (trace 
of the) magnetic field squared. 

There are also derivative terms for the loops them- 
selves: 

JC 2 (L) = 1 V Cs.s'.s" Re (dds) (dds') i S » , 

y S,S',S"eTl 

(19) 

where the representations are constrained so that the to- 
tal Z(N) charge of each term vanishes, es+es'+es" = 0, 
modulo N. 

Lastly, there are terms involving two derivatives of L, 
coupled to a field with Z(N) charge minus two: 

fC 3 (L) = ^ tr Re (AL) 2 Cs is , e 5 = -2 . 

5 Sell 

(20) 

For three colors, this series starts with the triplet loop. 

In the deconfined phase, the free energy is of order 
~ N 2 . Since all loops are of order one at large N, in the 
action we always assume that the potential V is multi- 
plied by an overall factor of ~ jV 2 at large N. The mass 
m 2 is then naturally of order one at large N; how the cou- 
pling constants Ks scale with N is discussed in Sec. IIII El 
For the kinetic terms, to contribute ~ iV 2 , the coupling g 
should scale like the gauge coupling, with g 2 N and g 2 N 
fixed at large N; then the couplings Cs,S',S" and Cs are 
both - 1. 

The renormalization of the loop potential was dis- 
cussed in The Wilson line in a given, irreducible 
representation undergoes mass renormalization, m^"; on 
a lattice, this is a power series in the coupling constant 
times the lattice spacing, 1/a. For a Wilson line at a 
temperature T, the renormalization constant, Z-r, and 
the renormalized loop, tfi, are given by 

in = in/Zn , Zn = exp(-m%°/T) . (21) 

Loops in irreducible representations do not mix under 
renormalization. 

The renormalization constants for the kinetic terms 
are similar, but more involved, than those for the po- 
tential. All kinetic terms undergo mass renormaliza- 
tion, with renormalization constants which are the ex- 
ponential of a divergent mass times the length of the 
loop. For loops without cusps, there is no condition to 
fix the value of these renormalization constants at some 
scale, while loops with cusps do require such a condition 
IE 0- F° r kinetic terms, additional renormalization 
constants, and conditions to fix their value at some scale, 
may be required. For example, each term in l|18|) experi- 
ences mass renormalization. In addition, since the term 
~ 1/g 2 arises from the electric field in the bare action, 
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at least in four spacetime dimensions, it will require an 
additional renormalization constant related to coupling 
constant renormalization. 

We have written down all possible terms involving two 
derivatives, but this classification may well be overly 
complete. In particular, renormalization will greatly re- 
strict the possible terms. The entire list is relevant in 
two spatial dimensions; the associated (3- functions can be 
computed in the ultraviolet limit |49| . In three spatial 
dimensions, a scalar field has dimensions of the square 
root of mass, and most of the above terms are non- 
renormalizable, and so can be ignored. 

Terms which arise from the coupling to fields which 
are not Z(N) invariant, such as quark fields, can also 
be included. The simplest possible coupling involves the 
trace of the chiral field: 



tr U Rc £ 



N 



(22) 



This, however, is chirally suppressed, proportional to m 2 . 
There are couplings which are not chirally suppressed, 
but these necessarily involve derivatives of the chiral field; 
these start as 



tr \d^U\ 2 Re. 



N 



(23) 



We do not know if this chiral suppression significantly 
affects the breaking of Z{N) symmetry; in mean field 
theory, this can be analyzed by using a nonlinear sigma 
model for the chiral fields [5(|. For an analysis of 
Z(N) breaking terms in terms of linear sigma models, 
see Mosey, Sannino, and Tuominen |3l| . 



III. INFINITE N 

A. Gross— Witten Point 

The effective action for loops is much more involved 
than for sigma models. Because loops have a potential, 
however, we can systematically perform a large N expan- 
sion by minimizing an effective potential. Corrections 
in l/N arise from fluctuations, which arise when spatial 
derivatives are included. In the rest of the paper, we ig- 
nore fluctuations, and the kinetic terms of Sec. Ill CI to 
concentrate on what is a matrix- valued mean field theory. 
Sec. Ill treats this matrix model at infinite N; Sec. IV, 
finite N. Even at infinite N, we stress that because we 
cannot compute the loop potential from first principles, 
all we are doing is characterizing all possible transitions 
in terms of couplings in the loop potential, the Ks of 1)170. 

If fluctuations are neglected, we can assume that the 
Wilson line is the same at each point in space, L(af) = L. 
Thus the functional integral reduces to a single integral, 



dL exp (-A 2 V(L)) 



(24) 



At large N, factorization implies that for arbitrary rep- 
resentations, normalized loops are just products of the 



fundamental, and anti- fundamental, loop [3j. Then the 
loop potential reduces simply to a power series in |^at| 2 , 



V(L) = -m 2 \£ N \ 2 + K4 (\£n\ 2 ) 2 + k 6 (\£n\ 2 ) 3 + - 



(25) 



We have relabeled the couplings ns in ltT7|l as K 2n , where 
the subscript now denotes the power of In- The U(l) 
symmetry is broken to Z(N) by a term ~ kjv Re (£n) N , 
but this is negligible at large N. 

As the potential is multiplied by an overall factor of 
~ A 2 , and the loops are normalized to be of order one, 
the natural guess is also to take the couplings Kg of or- 
der one at large N. 



As we discuss in Sec. IIII El this 
choice of couplings is far from innocuous. The advantage 
is that this assumption vastly simplifies the analysis, and 
allows us to gain insight which is otherwise obscured by 
technical complications Further, as far as the bulk 
thermodynamics is concerned, the only novel behavior 
emerges in this regime; otherwise, the deconfining tran- 
sition is rather ordinary. 

To find the stationary points of the functional integral, 
we introduce a delta-function into the integral [3^ as 



dL 



dX S(X - £ N ) exp (-A 2 V(A)) . (26) 



Here A is a (complex) number, equal to the value of £n = 
tr L/N for a given matrix L. As such, in the action 
we can replace V{£n) by V(A). The constraint is then 
exponentiated by introducing a field oJ, 



dX duJ dL exp (-A 2 V C 



,. = V(A) + <57(A - £ N ) 



(27) 



In this expression, both A and To are complex numbers, 
not matrices, or traces thereof. Usually, iTD is real at a 
stationary point, and so we define u> = itJ. 

The full integral is over A, 57, and L. We break off the 
integral over L, and define that piece as 



Zgw{u) 



J dL exp (A 2 lo Rc 



£n) 



(28) 



By an overall Z{N) rotation, we can choose any conden- 
sate for £pf to be real and positive. Both A and u> are 
complex fields, but to look for a real stationary point in 
£n, we need only consider the real parts. 

The integral Zqw was done by Brezin et al. [38| , Gross 
and Witten [3^ , and Aharony et al. . For the unitary 
matrix L, all that matters arc the eigenvalues of L, L^j = 
Sij exp(iQ!j). In the limit of infinite A, the number of 
eigenvalues is infinite, and it is convenient to introduce 
the density of eigenvalues, p(a). The solution for this 
density is: 



p(a) 



1 

2^ 



(1 + u) cos a) 



to < 1 ; 



(29) 
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p(a) 



1 at 2 a\ !/ 2 

— cos — 1 — uj sin 
7T 2 V 2> 



w > 1 



(30) 



where the latter only holds for uj sin 2 (a/2) < 1. 

We now make the following observation. The value of 
the condensate is 



da p(a) cos a . 



(31) 



With these expressions for the eigenvalue density, it is 
not difficult to show that 



and 



4>M 



£ (u) = 1 - 



uj < 1 



2^ 



uj > 1 



(32) 



(33) 



To implement these relations between the expectation 
value and to, we introduce yet another constraint field, £, 
along with a function, Vvdm(£)> as follows: 



Z G w(u) = / d£ exp (N 2 {uj£ - Vvdm{i))) 



(34) 



I is just a complex number, so the stationary point of 
this integral occurs when 



dV 



to 



Vdm 



or 



It is then trivial to see that if we choose 

f 

2 ; 



V V dm{t) = + r , £< 



(35) 



(36) 



V Vdm (l) = - l -\og{2{l-£))+ 1 l , l>\. (37) 

then l|35(l gives the desired relations, l|32|l and i|33|l . One 
can also check that (|34f) gives the correct result for 

Using this form of Vvdm iu Q34|). and replacing u = iu>, 
the original partition function becomes 

Z = I dX I duJ I dt exp (-N 2 V ) , 



V' = V{X) + iuJ(\-£) + Vvdm(£) ■ (38) 
Doing the integral over ZU fixes X = £, and leaves 

Z = J dl exp(-N 2 V eff {£)) . (39) 

The effective potential, 

V eff (£) = V(£) + V V dm{£) , (40) 



is the sum of the loop potential, which we started with, 
and Vvdm- Unlike the original integral over the matrix 
L, (|3*3|) is just an integral over a single degree of freedom, 
£. Since the effective potential is multiplied by an overall 
factor of TV 2 , at large N the true vacua of the theory are 
the stationary points of V e ff (£) ■ 

When we introduce Vvdm(£) m 03}, mathematically 
this is just the Legendre transformation of Zqw |35| . 
treating u> as an external source for £. This is nontriv- 
ial only because in Zqw, the measure for the matrix L 
includes the Vandermonde determinant. For this reason, 
we refer to Vvdm{£) as the "Vandermonde potential". 

Note that the Vandermonde potential is just the Leg- 
endre transformation of Zqw, and not of the full inte- 
gral. One could compute the Legendre transformation 
of Z, but of necessity, this would also include the loop 
potential. When the transition occurs, however, the ef- 
fective potential has degenerate vacua, and V e ff {£) is not 
a monotonic function of £. This leads to well known ambi- 
guities in the Legendre transformation |35|. In contrast, 
the Vandermonde potential is a monotonically increas- 
ing function of £, and there is never any ambiguity in its 
Legendre transformation. This remains true at finite N, 

Sec. nsu 

The Vandermonde potential vanishes at £ — 0, and di- 
verges, logarithmically, as I — > f. At £ = |, the value 
of Vvdm(£), and its first and second derivatives, are con- 
tinuous. The third derivative is not, since it vanishes for 

l 



£ < i, and is nonzero for 



> 



For a lattice theory, to is an inverse coupling constant. 
With a Wilson action, there implies a transition, of third 
order, in uj. This disappears if the lattice action is other 
than Wilson, however. In contrast, the Vandermonde 
potential always has a discontinuity, of third order, at 
£ = h. Because it is a transition of such high order, 
it may not affect the bulk thermodynamics significantly, 
but it is always there. 

The form of the Vandermonde potential is of interest. 
While it is nonpolynomial when £ > 4, ~ log(I — £), 
this is just eigenvalue repulsion from the Vandermonde 
determinant. This is why there are no constraints on the 
signs of the couplings in the loop potential: eigenvalue 
repulsion never lets the value of £ exceed one. 

About the origin, the Vandermonde potential is just a 
mass term, but in fact that is the most interesting thing 
about it. All terms of order ~ £ 4 , ~ £ e , and so on, 
conspire to cancel [lit . In a miracle of group theory, the 
leading correction is of very high order, starting out as 
~ £ N . 

Given the effective potential, it is then immediate to 
read off the phase diagram of the theory. For a given 
value of m 2 and the coupling constants, one just varies 
with respect to £, as one would for any other potential. 

For a given value of m 2 , we denote the stable minimum 
as £q. The confined vacuum has £o — 0: if the confined 
vacuum is a stable minimum of the loop potential, then 
because of the simple form of the Vandermonde potential, 
it remains so for the effective potential, with V e //(0) = 0. 
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A vacuum in the deconfined phase satisfies 

dVeff(l) 



Of 



= 



(41) 



with £q =/= 0. The transition occurs when the deconfined 
phase is degenerate with the confined phase, 



V e //(^o + ) = 



(42) 



£q denotes the value of the order parameter at the tran- 
sition, approaching it in the deconfined phase. 

The simplest possible example is to ignore all couplings 
in the loop potential, and simply take 



v{i) = -to 2 £ 2 

About the origin, the potential is 



V eff (£) = (-m 2 + l)£ 2 



>4 



and 

Veff(t) 



-to 2 £ 2 - 



log (2(1 -*)) + - 



£ > - 
~ 2 



(43) 



(44) 



(45) 



In the confined phase the potential is just a mass term. 
Its sign shows that a transition occurs at m 2 = 1, be- 
tween a confined phase, for m 2 < 1, and a deconfined 
phase, for to 2 > 1. In the deconfined phase, the conden- 
sate is the solution of 



which is 



-2m 2 £ 



k = - 
2 



1 



2(1-4) 



= 



(46) 



(47) 



In the confined phase, 

m 2 // = 2(1 - m 2 ) , to 2 < 1 . (49) 
In the deconfined phase, using (|46(l we can write 



_ 24-1 1 
e// " 24(1 - to) 2 ' 2 



About the transition, 



L eff 



4^ 



m- 



1 



r 



(50) 



(51) 



This difference occurs because at the Gross-Witten 
point, the new minimum is right at the point where there 
is a discontinuity of third order. Both masses vanish 
at the transition, but do so asymetrically, with different 
powers of |to 2 — 1|. 

To go further, we need to make some assumption about 
the relationship between the coefficients of the loop po- 
tential and the temperature. We assume a mean field 
relation between the adjoint mass and the temperature 
to 2 — 1 ~ T — Td, neglecting the variation of the cou- 
pling constants with temperature. As the transition is 
approached in the confined phase, m 2 ff ~ Td — T as T — > 
, while in the deconfined phase, m 2 e ^ ~ (T — Td) 1 / 2 
for T^T+@. 

At the transition, the potential vanishes at both de- 
generate minima, V e //(0) = V e ff(^o ) = 0. The derivative 
with respect to m 2 is discontinuous, though, vanishing in 
the confined phase, but nonzero in the deconfined phase, 



3V eff {£) 



dT 



dV eff (£) 



dm 2 



1 

4 ' 



(52) 



(There is another solution to the quadratic equation, but 
it has £q < i, and so doesn't matter.) 

The effective potential at the transition, to 2 = 1, is 
shown in Fig. ^ We call this the "Gross-Witten point" . 
At the transition, the potential vanishes for £ between 
and £q = |; it then increases for £ > i. Away from the 
Gross-Witten point, the potential is rather ordinary. In 
the confined phase, it is monotonically increasing from 
I = 0. In the deconfined phase, the global minimum has 
4 ^ 0. (For to 2 > 1, £ = is actually metastable, since 
although there is no barrier, the first derivative of V e // 
vanishes at I = 0.) As to 2 increases in value from one, 
the minimum moves from i to larger values, approaching 
one only asymptotically in the limit of m 2 — > oo. For any 
value of to 2 , either positive or negative, there is only one 
stable minimum. 

We can also compute the effective mass squared. With 
this method, it is simply the second derivative of the 
potential about the minimum, 



l eff 



d 2 Ve ff (e) 

d£ 2 



(48) 



which shows that the latent heat is nonzero, arising en- 
tirely from the deconfined phase. The transition is "crit- 
ical" first order [sj: although masses vanish, the order 
parameter jumps. 

Our discussion of the Vandermonde potential follows 
that of Aharony et al. ^3 • They argue that if the loop 
potential contains operators such as tr L 2 /A, V e f/(£) can 
only be plotted for I < i. For potentials which are sim- 
ply powers of the fundamental loop, however, we assert 
that the effective potential is meaningful, and most illu- 
minating to plot, over its entire range, for I : — > 1. In 
Sec. IIII El we outline how to compute the effective poten- 
tial for arbitrary loop potentials. 



B. Nonzero Background Field 

Minimizing the effective potential is an elementary 
exercise in algebra. We begin with a background 
Z(N) field, taking the loop potential to be 



V(L) 



-h Re i 



N 



TO 2 \£ 



(53) 
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In this case, there is a nontrivial minimum of the effec- 
tive potential for any nonzero value of ft. When £ Q < |, 
it is 



4 = 



2(1 -m 2 ) 



(54) 



while for Iq > 4, 




(55) 



Consider the case where h is infinitesimally small. If 
m 2 — 1 is not ~ ft, we can expand in ft: £o(h) « ft/2 when 
£ a < \ - 0(h), while 



e (h)^e (h = o)+ -i + 



ft 

4to 2 



(56) 



for €o > \ + 0(h). Thus away from £o ~ ^! on both sides 
of the transition the shift in is small, of order ~ ft 
(see e.g. Fig. 2 in |2£|). When m 2 — 1 ~ ft, though, the 
solution differs by a large amount from that for ft = 0. 
There is a third order transition when £o equals ^. The 
value at which this happens is easiest to read off from 
1)54(1. This occurs for m 2 = 1 — ft, at which point the 
effective mass squared rn 2 jj ~ 2ft. 

For an ordinary first order transition, as the value of 
a background field increases from zero, there is still a 
nonzero jump in the order parameter. The jump only dis- 
appears for some nonzero value of the background field, 
which is a critical endpoint. For larger values, there is 
no jump in the order parameter, and masses are always 
nonzero. For the loop potential of i|53|) . however, any 
nonzero background field, no matter how small, washes 
out the first order transition. For any ft ^ 0, there is 
still a transition of third order, at which the masses are 
nonzero. 

These results give an elementary interpretation of the 
calculations of Schnitzer |2£}. He computed the parti- 
tion function for quarks in the fundamental represen- 
tation coupled to SU (oo) gauge fields, taking space to 
be a very small sphere. Fields in the fundamental rep- 
resentation break the Z(N) symmetry, and act like a 
background field ft ~ Nf/N, where Nf is the number of 
flavors. When m 2 < 1 - ft, by {3J! £ ~ ft ~ Nf/N; sub- 
stituting this back into 1531) . the potential, and so the free 
energy, is of order V e ff(£o) ~ h 2 ~ (Nf/N) 2 . The first 
order transition disappears for any ft ^ |42| . leaving 
just a transition, of third order, when £q passes through 



i|43 



C. Mean Field Approximation 



The analysis in Sees. ITiTaI and iHIBl is most trans- 
parent for understan ding the physics. Following Kogut, 
Snow, and Stone [3. l4ll l42l|. we show how to compute 



the partition function in another way. This was origi- 
nally derived as a mean field approximation, although 
for loop models, we show that it is equivalent to a large 
N approximation. We consider the case where the loop 
potential includes just an adjoint loop; we also add a 
background field, taking the loop potential to be that of 
(|53|1 . The case of more general potentials follows directly, 
and is addressed in Sec. IV for N — 3. 

In the expression of J37J|, A only appears in the loop 
potential and in the constraint. Thus before doing the 
matrix integral, we can extremize with respect to A. With 
the V of lO, this fixes A = (lu - ft)/(2m 2 ). We then 
define a potential directly from the matrix integral: 



(57) 



exp(-N z V(uj)) = dLexp(N z ujRe£ N ) 



This leaves one remaining integral, over lu. In fact, since 
the relationship between A and lu is linear, we can sub- 
stitute one for the other, and so obtain the mean field 
potential, 

V mf (£) = m 2 £ 2 + V(2m 2 £ + ft) . (58) 

To follow previous notation, we have also relabeled A as 
I. 

To understand the relationship to a mean field approx- 
imation, consider a lattice theory in which there are fun- 
damental loops on each site, coupled to nearest neighbors 



with strength 



The number of nearest neighbors 



also enters, but for our purposes, scales out. Start with a 
given site, and assume that on adjacent sites, all matrices 
have some average value, £ = tr L/JV. The action is then 
the coupling, m 2 , times the average value I, times the 
value for that site, or ~ 2m 2 £Retr L/7V altogether. In 
a background field, one also adds ~ ftRetr h/N. With 
lu = 2m 2 £+h, this is the integral of l|57() . and gives V. The 
remaining part of the mean field potential, m 2 £ 2 , arises 
by imposing the mean field condition that the expecta- 
tion value of Re L equal the presumed average value, £ 

Q 

Several aspects of the mean field approximation are 
much clearer when viewed as the large ./V expansion of 
a loop potential. We see that the mean field action, 
involving the real part of the fundamental loop, actu- 
ally arises from a loop potential with an adjoint loop, 
= to 2 |^jv| 2 . Thus if we add an additional term propor- 
tional to an adjoint loop, m 2 |^jv| 2 , to the mean field ac- 
tion, = 2m 2 £Ke£N, this is equivalent to a shift in the 
coupling for the adjoint loop, m 2 — > m 2 + rh? . This was 
proven previously in 3j, by more indirect means. 

For zero background field, the mean field potential is 



V m f(i) =m 2 {l-m 2 )£ 2 , £< 



2m 2 ' 



(59) 



and 



V m f (I) = -2m 2 £ + m 2 £ 2 + ^ log(2m 2 £) + | , I > ^-L . 

?60) 
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Notice that the mean field potential, V m f(i), does not 
agree with the previous form of the effective potential, 
V e ff{£) in lfHjl an d BB - Even the point at which the 
mean field potential is nonanalytic, £ = l/(2m 2 ), is not 
the same. 

Although the mean field and effective potentials differ, 
they do give the same vacua. In the deconfined phase, 
the stationary point of H60[l satisfies 



-2m 2 (l-^)+2^ = 0, 



(61) 



which agrees with (|46|l . One can also check that for h ^ 
0, the solution coincides with l|54|) and (|55|l . 

Second derivatives of the mean field and effective po- 
tentials are not the same, even at the correct vacuum £$. 
For instance, in the confined phase the second derivative 
of {23} is = 2m 2 (l - to 2 ), not = 2(1 - to 2 ) of gU). They 
agree to leading order in 1 — to 2 about to 2 — > 1, but not 
otherwise. Likewise, in the deconfined phase, the second 
derivative of V m f only agrees with that of V e //, {ST]) , to 
leading order in y/m 2 — 1 as to 2 — > 1. 

To understand this discrepancy, we define a mass not 
by the potential per se, but by the response to a back- 
ground field. Computing the partition function in the 
presence of h ^ 0, the second derivative is 



1 d 2 
Is^dh 2 



log Z{h) 



((Re £ N ) 2 ) - (Re £ N ) 2 . (62) 



h=0 



The great advantage of the approach of Sec. IIII Al is 
that the background field h only appears linearly in the 
action, through the loop potential. In this case, l|62|) = 
— 1/(2A^ 2 to 2 j^), where rn 2 ^ is just the second derivative 
of the effective potential, (J2SJ. 

In the mean field approach, however, h enters into 
V(2to 2 ^ + h); this function has terms of both linear and 
quadratic order in h. In order to compute l]62p. then, it 
is necessary to include the terms quadratic in h. This 
can be done, but is tedious. The result is that the effec- 
tive mass, defined properly from H62[) and computed with 
V m /, agrees with that obtained so easily from V e //, f4l?|) 
and ijSTJjl . 

In the presence of a nonzero background field, following 
Damgaard and Patkos [43j we have checked that there 
is a transition of third order when £ passes through \ . 
Even this is not obvious for the mean field potential, since 
the point at which V m f is nonanalytic depends upon the 
value of to 2 . 

In the end, the two methods must agree for physical 
quantities. One is only doing the integrations in a differ- 
ent order, and there are no subtleties in the integrands. 
On the other hand, in the next section we see that at 
finite N, the mean field approach is advantageous for 
both formal calculations in the confined phase, and for 
numerical computation. 



D. Away from the Gross— Witten Point 

We now return to the case of zero magnetic field, and 
consider higher interactions. Consider a quartic potential 

El, 



V(L) 



N 



Ki (\£ 



(63) 



It is immediate to compute the phase diagram. There 
is a line of second order transitions for K4 > 0, along 
—to 2 + 1 = 0, and a line of first order transitions for 
K4 < 0, along some line where —to 2 + 1 > 0. These 
meet at what appears to be a tri-critical point, where 



K4 



1 = 0: this is the Gross- Witten point. 



It is not a typical tri-critical point, however. If this 
were an ordinary scalar field, then £q , the value of the 
order parameter just above the transition, would be zero 
along the second order line, zero at the tri-critical point, 
and then increase from zero as one moves away from the 
tri-critical point along the first order line. 

For the matrix model, however, £q = along the sec- 
ond order line, but then jumps — discontinuously — to 
^0 = I at tne Gross- Witten point. 

For example, for small and negative values of K4, the 
transition occurs when 



1 



l«4| 

4 



At the transition, £ q increases from i as 



(64) 



(65) 



The masses are always nonzero; at the transition, in the 
confined phase, 



m lff 



l«4| 

2 



while in the deconfined phase, 



(66) 



(67) 



When K4 — > 0, the mass in the confined phase vanishes 
more quickly than in the deconfined phase. This is simi- 
lar to what happens when K4 — 0, (|49|l and l|51l) . 

Our analysis differs from that of Aharony et al. [l8| . 
The caption of their Fig. 3 is correct, stating that for 
small K4, the transition occurs at l|64|) . In the Figure, 
however, the curve for this to 2 is actually in the decon- 
fined phase. This is not evident, since they only plot 
the potential for £ < i. It can be seen by noting that 

at the transition, V e ff(£~o ) = 0; as £q > \, il6"5|l . then 
Veff(h) > 0, and not zero, as in their Figure. 

Moving to increasingly negative values of K4, the value 
of £q increases as well, as do the masses. For example, 
consider K4 = —1. Numerically, we find that the transi- 
tion occurs when to 2 w 0.456121..., with £q w .841176.... 
The effective potential at the transition is illustrated in 
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Fig. [5] The masses are very different in the two phases: 
in the confined phase, ni^j = 1.08776..., while in the de- 
confined phase, m 2 e ff = 10.4184.... Thus, for such a large 
value of — K4, the transition appears to be a perfectly or- 
dinary, strongly first order transition. The masses are al- 
ways nonzero, including at the transition point. In Fig-El 
the potential does have a discontinuity, of third order, at 
I = |; however, as £q > \, this is of little consequence. 

We remark that because the masses are very different 
in the two phases, the potential in Fig. |3 cannot be ap- 
proximated by a quartic potential in £. We return to this 
for N = 3 in Sec. HVBl 

As K4 — ► — oo, m 2 w |k 4 |, and £q w 1 — 1/4| K4I . Due 
to eigenvalue repulsion in the Vandermonde potential, 
the expectation value is always less than unity. Since 
£^ — ► 1 as K4. — ► — co, in this case the deconfined phase 
is arbitrarily close to a truly perturbative gluon plasma 
from temperatures of on up, with the transition as 
strongly first order as possible. It is reasonable that what 
is a strong coupling phase in the loop model corresponds 
to a weakly coupled regime in the underlying gauge the- 
ory. 

Next, consider adding six-point interactions, ~ 
«6(K-iv| 2 ) 3 - F° r the purposes of discussion, we take k 6 to 
be positive, although this is not necessary. For Kg > 0, 
the phase diagram, in the plane of m 2 and K4, is greatly 
altered. There is a second order line when K4 > 0, which 
meets a first order line for K4 < 0, but they meet at a true 
tri-critical point, for K4 = 0. At this tri-critical point, the 
jump in the order parameter vanishes, as does the mass. 

With only quartic interactions, either £q — or 1$ > | 
|18|. This is no longer true when six-point interactions 
are included. Consider the loop potential 

V(L) - k 6 \£ n \ 2 (\£ N \ 2 - i 2 c ) 2 + (m 2 - l)\£ N \ 2 , (68) 

where £ c is some number. If £ c < i, then there is a first 
order transition when m 2 = 0, as £q = £ c at the tran- 
sition. If £ c > 5, then one must use the Vandermonde 
potential in (|37[1 . and while the transition is still of first 
order, it does not occur when m 2 — 0. We do not work 
out specific examples, since we only wanted to make the 
point that it is possible to have first order transitions in 
which £q assumes any value between zero and one. 

In the cases in which there is a second order transition, 
such as K4 > 0, there is nothing remarkable about it: it 
is in the universality class of a U(l) spin, and exhibits 
standard critical behavior. 

In all cases, if £q < i, then there is a third order tran- 
sition when the value of £0 passes through i 18]. Since 
the masses are nonzero at the transition, this third order 
transition does not appear to be of much consequence. If 
£q > 5, there is no third order transition, just one first 
order transition. 

The only way to have a first order transition, in which 
both masses vanish, is if all couplings vanish: K4 — — 
... = 0. For this reason, we term the Gross- Witten 
point an ultra-critical point, since an infinite number of 



couplings must be tuned to zero in order to reach it. 

E. Couplings at Large N 

In this subsection we comment on how couplings in the 
loop potential scale with N. 

To one loop order in the perturbative regime, terms 
such as as terms ~ tr L 2 /jV arise in the loop poten- 
tial [H E E |H HI El EH This is ~ l at large 
N, and so allowed. In terms of normalized loops, how- 
ever, tr L 2 /N ~ N(£ N 2 — £ 2 N ), where i N 2 is a loop 
for a two-index tensor representation, Eq. (42) of 
Thus |tr 1?/N\ 2 ~ N 2 \£ m - l 2 N \ 2 . Traces such as 
I tr L p /iV| 2 , which are also allowed in the loop potential, 
are ~ 7V 2 ( p_1 ) times differences of normalized loops. 

For the general loop potential of 117fl . the Gross- 
Witten point occurs when m 2 = 1, with all other cou- 
plings ks = 0. Thus we studied a weak coupling regime 
close to the Gross- Witten point, where the ks ~ 1 at 
large N. A regime in which the couplings Kg grow with 
powers of N is one of strong coupling, far from the Gross- 
Witten point. 

For strong coupling, the loop potential depends not 
just upon the fundamental loop, but on all N — 1 de- 
grees of freedom of the matrix L. A convenient choice 
would be the fundamental loop, tr L/iV, plus tr L 2 /A, 
tr L 3 /N, . . . tr L^" 1 /N . It is then necessary to compute 
all potentials in the associated N — 1 dimensional space. 
This will be involved, since an expectation value for the 
fundamental loop automatically induces one for tr L 2 /A, 
etc. 

Aharony et al. studied this strong coupling regime by 
adding |tr L 2 /A| 2 to the loop potential |18|. Using the 
solution of Jurkiewicz and Zalewski |40| . they find that 
the behavior is similar to that in weak coupling when 
K4 7^ 0. The bulk thermodynamics only exhibits ordi- 
nary first or second order transitions, plus third order 
transitions when £q passes through | . This is most likely 
valid everywhere in the strong coupling regime. 

Thus whether or not deconfinement is near the Gross- 
Witten point when N > 3, in the space of couplings Kg, 
it is manifestly the most interesting place it could be. 

IV. FINITE N 
A. Z(N) Neutral Loops in the Confined Phase 

What made our treatment at infinite N so simple is the 
assumption that because of factorization, the only (nor- 
malized) loops which enter are those for the fundamental 
representation. At finite N, one cannot avoid considering 
loops in higher representations. 

We consider only effects of the loop potential, neglect- 
ing kinetic terms. While this is exact at infinite N, it 
is not justified at finite N. Thus the following anal- 
ysis should be considered as the first step to a com- 
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plete, renormalized theory, including the kinetic terms 
of Sec. HTT1 

In this subsection we begin with some general observa- 
tions about the expectation values of loops in arbitrary 
representations: 

(£ K ) = JdL £ n exp(-A 2 V(L))/Z , (69) 

with Z that of (|2~4"|l . 

We concentrate on the stationary point of the partition 
function. This classical approximation is exact at infinite 
N; how well it applies at finite N, even to constant fields, 
is not evident. For N = 3, however, the overall factor of 
N 2 [or in fact, N 2 — 1 = 8, cf. eq. (|7S|l ] in the exponen- 
tial suggests that this might be reasonable. While less 
obvious for N — 2, we suggest later a very specific test 
which can be done through numerical simulations. 

In fact, we are really forced to consider just the sta- 
tionary point of the partition function. If we did the com- 
plete integral, then at finite N the Z(N) symmetry would 
never break spontaneously. In a phase which we thought 
was deconfined, all Z{N) transforms of a given vacuum 
would contribute with equal weight, so that in the end, 
all Z{N) charged loops would vanish (Z(N) neutral loops 
would be nonzero). This is just the well known fact that a 
symmetry only breaks in an infinite volume, or at infinite 
N ^^|. By taking the stationary point of the integral, 
we are forcing it to choose one of the vacua in which the 
Z{N) symmetry is broken. 

As noted before, the matrix L depends upon N — 1 
independent eigenvalues. Instead of the eigenvalues, we 
can choose the N — 1 independent quantities to be tr L, 
tr L 2 , and so on, up to tr L^ -1 . 

We prefer to deal with a set of loops. The first 
trace is of course the fundamental loop, £n = tr L/iV. 
Any potential is real, and so we also choose the anti- 
fundamental loop, £jf = l* N . As a Young tableaux, the 
anti-fundamental representation corresponds to the anti- 
symmetric tensor representation with N—l fundamental 
indices. 

There is no unique choice for the remaining loops, but 
as the anti-fundamental representation is a purely anti- 
symmetric representation, we can chose to work with only 
anti-symmetric representations. Thus we choose loops 
for anti-symmetric tensor representations with j indices, 
j = 1 to N — 1, denoting this set as TZa- 

For N = 2, TZa is just the doublet loop, which is real. 
For TV = 3, it is the triplet and anti-triplet loops. For 
-/V = 4, TZa includes the quartet, anti-quartet, and sextet 
loops. The quartet transforms under a global symmetry 
of Z(4), while the sextet only transforms under Z(2); 
both are needed, to represent the case in which only Z(2) 
breaks, but Z(A) doesn't. Including all of the loops in 
TZa ensures that all possible patterns of Z(N) symmetry 
breaking can be represented for N > 4. 

The previous form of the loop potential, l|17fl . is an 
infinite sum of loops in Z(N) neutral representations, the 
TZq. Instead, we now take the potential to be a function 



only of the loops in TZa, 

V(L)=V(£ S ) , SeTZ A . (70) 

Previously, the potential was a sum over the loops in 
TZq, each appearing only to linear order. Now, only the 
N—l loops in TZa enter, but these do so through all 
Z(N) invariant polynomials, to arbitrarily high order. 

This is useful in introducing constraints into the par- 
tition function. One might have thought that it is neces- 
sary to introduce constraints for all fields which condense, 
but in fact, it is only necessary to introduce constraints 
for fields which appear in the action |35J. Even so, if the 
action is written in terms of the TZq, there are an infinite 
number of such loops. The advantage of the TZa is that 
we only have to introduce N—l constraints. 

As usual, we start by introducing constraints into the 
partition function, 

Z= f dL H fdXs S(X s -£s) exp(-A^ 2 V(A 5 )) . 

J SeTZ A J 

where in the potential we have used the delta-functions 
to replace V(£s) by V(As). Introducing constraint fields, 

Z= j dL Y[ fd\ s I duj s cxp (-N 2 V cons .) , 

J Sell A 

(72) 

we obtain the constraint potential, 

V cons . = V(X S ) + £ iusiXs- £s) ■ (73) 

(In this subsection alone we write u>s instead of Us, since 
all ojs's vanish at the confined stationary point.) We now 
follow the mean field approach, defining the potential 

exp{-N 2 V(uj s )) = j dL expiN 2 £ u; s £ s ) . (74) 
After integrating over L, we obtain 




(75) 

where 

V m /(A 5 ,cJs) =V(A 5 )+ iu s X s + V(iuj s ) . (76) 

SeTZ A 

This is the mean field potential, written in terms of the 
TZa loops. 

By an overall Z(N) rotation, we can always assume 
that the stationary point is real. Thus in order to deter- 
mine just the stationary point, only the real part of any 
loop will enter. The integral which determines V for the 
fundamental loop is familiar : see Table 12 of [3|| and 
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|47|. The integral in (|74|l is more general, involving all 
loops in TZ_a- 

We use this formalism to prove an elementary theo- 
rem which is mathematically trivial, but physically im- 
portant. The obvious guess for the confined phase is 
where the expectation values of all loops in the action 
vanish, As = for all S in TZ_a. We can always de- 
fine the loop potential so that it vanishes when all £s 
vanish; any constant term in V drops out of the ratio 
in 1)69(1. Further, as the potential is Z{N) neutral, and 
as all \ s carry Z(N) charge, the first derivative of V 
with respect to any As has Z(N) charge, and vanishes 
if As = 0. Thus all constraint fields vanish at the sta- 
tionary point, lus = 0. Z(N) invariance also implies 
that the first derivative of V(u>s) with respect to any ws 
vanishes. Altogether, V m f = in the confined phase, 
As = us = 0, and this point is extremal. If the poten- 
tial vanishes, though, the expectation value of any loop 
is simply an integral over the invariant group measure. 
For any (nontrivial) representation, however, whatever 
the Z(N) charge of the loop, its integral over the invari- 
ant measure vanishes identically, J £fi dL = 0. Conse- 
quently, the expectation values of all loops vanish: 

(in) = , T<T d - (77) 

While the confined vacuum is extremal, stability is de- 
termined by second derivatives of the potential, and only 
holds for T < TJ. 

This is very unlike what would be expected merely on 
the basis of Z(N) symmetry. For Z(N) neutral loops, 
such as the adjoint, there is no symmetry which pro- 
hibits their acquiring a nonzero expectation value in the 
confined phase. We stress that this result is valid only in 
mean field theory, when fluctuations from kinetic terms 
are neglected. 

For three colors, lattice data indicates that the expec- 
tation value of the renormalized adjoint loop is very small 
below Td Q • Since the deconfining transition is of first 
order for three colors (in four spacetime dimensions), this 
is not the best place to look for the expectation value of 
Z(N) neutral loops below Td- Rather, it is preferable to 
study a deconfining transition of second order, as occurs 
for N = 2 in three or four spacetime dimensions, and 
even for N = 3 in three spacetime dimensions. We re- 
mark that for a second order transition, universality im- 
plies that two point functions, such as (l a d{x)\(-N(0)\ 2 ) , 
scale with the appropriate anomalous dimensions. Uni- 
versality, however, places no restriction on the one point 
functions of Z(N) neutral loops. Merely on the basis of 
Z(N) symmetry, one does not expect them to be small 
in the confined phase: this is a clear signal that the un- 
derlying dynamics is controlled by a matrix model, and 
not just by some type of effective Z(N) spin system. 

In this regard, Christensen and Damgaard, and 
Damgaard and Hasenbusch [44], also noted that in a 
classical approximation, not only do the expectation val- 
ues of Z{N) neutral loops vanish in the confined phase, 
but for a second order deconfined phase, they vanish 



like a power of the fundamental loop. As T — > T^, 
(f-n) ~ (£n) p+ (£*n) p ~ , up to corrections involving higher 
powers of the fundamental and anti-fundamental loops. 
Thepowers are identical to those of factorization at large 
N |3j, but the coefficient is not the same; to satisfy factor- 
ization, the coefficient is one, up to corrections of order 
l/N. For AT = 2, it is 2/3. 

We conclude this section with a comment. There is al- 
ways a given element of SU (N) for which the trace of the 
fundamental loop vanishes; e.g., in SU(2) it is the diag- 
onal matrix ±(1, —1). Thus one might imagine modeling 
the confined vacuum as an expansion about such a fixed 
element of the group group [33. Because these matrices 
are not elements of the center of the group, however, they 
are not invariant under local SU(N) rotations: they rep- 
resent not a confined, but a Higgs, phase. This is very 
different from the above, where in the confined phase, one 
integrates over all elements of the gauge group with equal 
weight. While the expectation value of the fundamental 
loop vanishes in a Higgs phase, those of higher represen- 
tations do not; that of the adjoint loop = — 1/(A 2 — 1). 
This does not agree with lattice simulations for N = 3, 
which find a very small value for the octet loop below Td, 
<l/8 0. 



B. Three Colors 

For N = 3, the analytic solution of large N is not avail- 
able. There are various approximation schemes which 
one can try. 

For example, one can ex pan d about the confined phase, 
I = 0. From Eq. (2.22) of [4l|. however, one can see that 
while such an expansion works well for small I < 0.2, 
it fails at larger values, and so is not of use near the 
Gross- Witten point. (Note that absorbs the number 
of nearest neighbors into their parameter a, so a = 6£; 
also, their = (N 2 - \)m 2 /N 2 = 8m 2 /9 due to our 
definition of the prefactor of V in eq. (JTSJ.) 

One can also expand about the deconfined phase, tak- 
ing to 2 — > 00, which forces £ — > 1. Since for N = 3 we 
need only consider the triplet loop, the integral for the 
mean field potential, (|74|) . is the N = 3 analogy of 15711 . 
At the stationary point, u> = 2m 2 £, so at large to 2 , where 
£ w 1, u> is also large. Mathematically, the integral for 
the mean field potential is the same as arises in strong 
coupling expansions of a lattice gauge theory. The result 
for large ui is given in Table 12 of [3g. The result starts 
as ~ lo, plus a term ~ log(w), and then a power series in 
l/ui. The first term is trivial, due to the fact that £ w 1. 
The term ~ log (a;) represents eigenvalue repulsion from 
the Vandermonde determinant. This was seen before at 
large N, as the term ~ log(l — £) in the Vandermonde po- 
tential, (|37|) . One finds that while this expansion works 
well near £ w 1, it does not appear useful for smaller 
values. 

The failure of these perturbative expansions is not par- 
ticularly remarkable. In terms of £, the Gross-Witten 
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point is identically midway between a confined, and a 
completely deconfined, state. There is no reason why an 
expansion about either limit should work, although they 
might have. Even so, for N — 3 the matrix model is just 
a two dimensional integral. The regions of integration 
are finite, and there are no singularities in any integral. 
Thus we can just do the integrals numerically. 
We define the partition function as 

Z = J dL cxp(-(iV 2 - l)V(L)) . (78) 

We multiply the loop potential by TV 2 — 1 , instead of N 2 . 
For the loop potential, this is a matter of convention, 
but for the Vandcrmonde potential, numerically we find 
that with this definition, the N = 3 results are closer to 
those of N — oo. This is not surprising: in expanding 
perturbatively about the deconfined state with t — 1, an 
overall factor of the number of generators, which for an 
SU(N) group is N 2 — 1, arises naturally 36]. 

By the previous section, we consider the loop potential 
as a function of the triplet and anti-triplet loops. The 
most general potential is then 

V(L) = -m 2 \e 3 \ 2 + « 3 Re (£ 3 ) 3 + k 4 (\£ 3 \ 2 ) 2 + ... (79) 

These terms represent octet, decuplet, and 27-plet rep- 
resentations of SU(3). Other terms of higher order in £ 3 
include one of pentic order, and two of hexatic order |'26| . 

We then mimic the analysis of Sec. IIII Al to obtain the 
Vandermonde potential, Vvdm(fy- The only difference 
is that at large TV, the expectation value of iff, in the 
presence of an external source, is known analytically, l|32|) 
and l|33|l . For N — 3, it is necessary to determine this 
relationship numerically, from the integral 

£{u>) = J dL Re £ N exp((N 2 - 1) w Re£ N ) , (80) 

where w is a real source. Given £{u>), it is then trivial to 
invert this function, to obtain the source as a function of 
the expectation value, w(£). Following l|35|l . the Vander- 
monde potential is then just the integral of the source, 
with respect to £: 

V V dm(£) = [ dt , (81) 

Jo 

As at infinite N, when N — 3 the Vandermonde po- 
tential is a monotonically increasing function of £: it van- 
ishes at I = 0, and diverges, logarithmically, as £ — > 1. 
Because this function is monotonically increasing, if we 
know £(ui), then there is no ambiguity in inverting it, to 
obtain uj(£). 

The numerical result for N — 3 is compared in Fig. 
to the analytical expression for N = oo from eqs. I)36|) 
and H37|) . The N — 3 result is always less than that 
for N = oo, but even up to £ « 0.8, they lie within 
a few percent of one another. This is most surprising: 
the U(l) symmetry of N — oo is broken at finite N to 



Z(N) by operators which start as £ N . For N — 3, this 
is an operator of low dimension, ~ £ 3 . Indeed, in Fig. El 
we also show the result for the N = 2 Vandermonde 
potential, and find that even that isn't so far from the 
N = oo result. 

To obtain the effective potential, we simply add the 
loop potential V to Vvdm- The result for the three-color 
analogue of the Gross-Witten point, i.e. the potential 
from eq. I|79(l with all couplings Ki = 0, is depicted in 
Fig.0|by the solid line. The value of the order parameter 
at the transition, £q « 0.485, is very close to the N — oo 
value of i |4l| . Since the Vandermonde potential for 
N = 3 is so close to N = oo, this small shift, by » 3%, 
is reasonable. Similarly, as can be seen graphically, the 
potential is almost flat at Td- Masses are small at the 
transition, and there is a small barrier between degener- 
ate minima of the potential. Physically, this corresponds 
to a nearly vanishing interface tension between phases. 

As mentioned in the Introduction, lattice data for 
SU(3) Yang-Mills indicates that 4(T+) is less than ±, 
perhaps w 0.4 ± 10% Q- Additional lattice simulations 
are needed to fix the jump in the triplet loop more pre- 
cisely. Nevertheless, in what follows we shall assume that 
£q = 0.4 to illustrate the effects of interactions. 

We consider in detail adding a cubic term to the loop 
potential, k 3 ^ 0. We have also considered adding higher 
terms, such as a quartic term, K4 7^ 0, but found numeri- 
cally that the results are rather similar. With k 3 w 0.146, 
we obtain £q w 0.4 for the triplet loop. The effective po- 
tential is shown in Fig. and exhibits an even smaller 
barrier between the two phases than for k 3 = 0. 

Numerically, we also find that at the transition, the 
masses in the two minima are equal, to within a few per- 
cent, for both K3 = and = 0.146. This allows us to make 
the following observation. Any polynomial approxima- 
tion to the potential fails for £ as 1, due to a logarithmic 
term which represents eigenvalue repulsion. However, we 
can always use a polynomial approximation for small £. 
Consider a polynomial in £ to quartic order, of the form 
~ £ 2 (£ — £q) 2 : this represents two degenerate minima, 
at £ — and £ = £$ . We find that at the transition, 
such a form is approximately valid for £ < 0.6. Further, 
it is clear that this potential gives equal masses in both 
the confined and deconfined phases, which is nearly true 
numerically. 

Such a quartic parametrization of the effective poten- 
tial is, in fact, the basis for the Polyakov loop model 

n .in in m m m m us eiii. Thevaiueof 

£q in the Polyakov loop model, ~ .55 |2J|, is close to that 
found from the lattice |3j,|j|. Moreover, the form of the 
potential is graphically very similar to that found in the 
matrix model, with an extre mely small barrier between 
the two phases: see Fig. 1 of |24|. 

We also computed the effective mass in the matrix 
model. At infinite N, the value about the Gross-Witten 
point is given in eq. I|50|l . For N = 3, the effective mass is 
obtained numerically by computing the curvature of V e // 
about the non-trivial minimum £q, cf. eq. I|48[). The re- 
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suits for K3 = and K3 7^ are shown in Fig. as a func- 
tion of the expectation value of the fundamental loop. At 
the Gross- Witten p oint £ p{T+) = \ and m 2 eff {T+) = 0, 
as discussed in Sec. IIII Al When TV = 3, at the transition 
the effective mass is small, but nonzero. In agreement 
with the discussion of V e // above, when K3 ^ 0, both the 
expectation value of the triplet loop, and the effective 
mass, are smaller than for K3 = 0. Except for this small 
effect, for both values of K3 the effective mass for N = 3 
is close to that at infinite N. This happens because the 
Vandermonde potential at N — 3 is close to N = 00, and 
K3 = 0.146 is a small value. For example, since £ w \ 
and to 2 sa 1 at the transition, the cubic coupling, k^£ 3 , 
is s» 8% relative to the adjoint loop, m 2 £ 2 . 

To compute the expectation values of loops in higher 
representations we use a mean field analysis, like that 
of Sec. rrncl Q. We checked that for n 3 = 0, the ex- 
pectation value of the triplet loop, as computed from 
the effective potential, agrees with the mean field result. 
The expectation values of the loops in the triplet, sextet, 
octet, and decuplet representations of SU(3>), computed 
in this way, are shown in Fig. [5] We plot them as func- 
tions of the coupling m 2 , divided by the critical value 
at which the transition occurs: m 2 rit w 0.91 for k 3 = 0, 
and m 2 rit w 0.97 for k 3 = 0.146. For k 3 ^ 0, the cubic 
term in the loop potential reduces the expectation val- 
ues of all loops, not just that of the triplet loop. When 
m 2 — > oo, the effect of the cubic interaction diminishes, 
as the expectation values of all loops approach that for 

Ka = 0. 

In agreement with the original results of Q], we find 
that the triplet loop, as measured on the lattice, agrees 
approximately with that for K3 = 0.146, assuming a lin- 
ear relation between m 2 and the temperature about 
The lattice data is not sufficiently precise to make a de- 
tailed comparison, however. In particular, as the loop 
approaches one, this is not an accurate way of determin- 
ing the relationship between m 2 and the temperature. 
We return to this point shortly. 

For representations beyond the triplet, the difference 
loop, introduced in [3j, is the remainder between the ex- 
pectation value of the loop, and the result expected in 
the large ./V limit, which is a product of fundamental 
(and anti- fundamental) loops: 

Sin = {in) . (82) 

The integers (p + ,p_) are (2,0) for the sextet, (1,1) for 
the octet, and (3, 0) for the decuplet representation. The 
difference loops are plotted in Fig. Here, too, we find 
that the cubic interaction only matters near the transi- 
tion: 5£g, S£q and S£±o are slightly smaller in magnitude 
near when K3 > 0, relative to their values for K3 = 0. 

In all cases, we find that the difference loops are much 
smaller than those measured on the lattice, Fig. 9 of 
To make a detailed comparison, it is necessary to know 
the precise relation between the parameter m 2 of the ma- 
trix model and the temperature. Even without this re- 
lation, however, the sextet difference loop is about four 



times larger on the lattice than in the matrix model. For 
the octet difference loop, the lattice data exhibits a sharp 
spike close to the transition, about six times larger than 
the value in the matrix model. There is, as of yet, no data 
for the decuplet loop from the lattice. To describe these 
behaviors in a matrix model, it will certainly be neces- 
sary to include fluctuations, due to kinetic terms in the 
action. It is not apparent whether the lattice data can 
be fit with couplings and kinetic terms whose coefficients 
are independent of temperature. 

We conclude by discussing the relation ship of our re- 
sults to the Polyakov loop model |H El El El El El 

ei ei in mm. In this model, the pressure is assumed 

to be a potential for £ times T 4 . As discussed above, this 
form does seem to work well near the transition. In order 
to fit the pressure away from T^, however, it is necessary 
to assume that the variation of m 2 is not simply linear 
in the temperature. As shown in p7|. by temperatures 
of « 3Td, the pressure is nearly a constant times T 4 ; this 
requires that to 2 is nearly constant with respect to tem- 
perature. As noted above, when the value of the loop is 
near one, and given the uncertainties in the lattice data, 
we cannot exclude such a change in the relationship be- 
tween m 2 and the temperature. 

For N = 3, assuming that the real part of the Polyakov 
loop condenses, one can compute masses for both the 
real and imaginary parts of the Polyakov loop |26j . The 
effective masses discussed above refer only to that for 
the real part of the Polyakov loop. This is a (trivial) 
limitation of the matrix model, and is easiest to un- 
derstand in the perturbative limit, when m 2 — > 00, so 
I — > 1. In this limit, L sa 1 + er, where a is an ele- 
ment of the Lie algebra, and the action is ~ m 2 tr a 2 . 
Integrating over a, by 162|) the effective mass is related 
to \ jm\ ss ~ J da (tr (a 2 )) 2 exp(— m 2 tr er 2 ), and gives 
m 2 jY ~ to 4 . This is the same behavior as at large N, 
as can be read off from (|50|l when £q — > 1. If we do the 
same for the imaginary part of the loop, which starts as 
Imag(tr L) ~ tr er 3 , we would conclude that its effective 
mass squared is ~ to 6 . This is wrong, however: in the 
full theory, the ratio between the effective mass for the 
imaginary, and real, parts of the fundamental loop must 
be 3/2, as it is in the gauge theory in the perturbative 
limit. The problem, as noted by Brezin et al. [38|, is sim- 
ply that because there are no kinetic terms, one cannot 
generally compute two point functions in a matrix model, 
but only one point functions. It is really exceptional that 
we could obtain even the mass for the real part of the fun- 
damental loop. Of course, this is automatically remedied 
by adding kinetic terms to the matrix model. 



V. CONCLUSIONS 

In this paper we developed a general approach to the 
deconfining transition, based upon a matrix model of 
Polyakov loops. The effective action for loops starts with 
a potential term. The Vandermonde determinant, which 
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appears in the measure of the matrix model, also con- 
tributes a potential term. 

At infinite N, the vacua are the stationary points of an 
effective potential, the sum of loop and Vandermonde po- 
tentials. Because of the Vandermonde potential, there is 
a transition with just a mass term, at the Gross- Witten 
point. This transition is of first order, but without an in- 
terface tension between the two phases. In the space of all 
possible potentials, which include arbitrary interactions 
of the Wilson line, the Gross- Witten point is exceptional. 
Away from it, the transitions appear ordinary, either of 
first order, with a nonzero jump of t and nonvanishing 
masses, or of second order, with vanishing masses and 
no jump in the order parameter. (There are also third 
order transitions, with no jump in the order parameter 
and nonzero masses.) Only by tuning an infinite number 
of couplings to vanish does one reach the Gross- Witten 
point, where the order parameter jumps, and yet masses 
are zero. 

We investigated the N = 3 matrix model about the 
Gross- Witten point, including interactions. We found 
that the N = 3 Vandermonde potential is very close to 
N = oo, so that for small interactions, the N — 3 effective 
potential strongly resembles that of infinite N. We found 
that in the N — 3 matrix model, corrections to the large 
N limit are small, ~ 1/N. 

What is so surprising about the lattice data is that 
the SU(3) deconfining transition is well described by a 
matrix model near the Gross- Witten point 0. In the 
confined phase, on the lattice the expectation value of 
the Z(3) neutral octet loop was too small to measure; 
in a matrix model, it vanishes identically. At the transi- 
tion, the renormalized triplet loop jumps to ~ 0.4 00, 
which is close to the Gross- Witten value of \. Further, 
masses associated with the triplet loop — especially the 
string tension — do decrease significantly near the tran- 
sition [Tol ]. 

Of course, the proximity of the SU(3) lattice data to 
the Gross-Witten point could be serendipitous, due to 
N = 3 being close to the second order transition of 
N = 2. For N > 4, the lattice does find a first order 
deconfining transition 0, 0, . The lattice data also 
shows that at a fixed value of T/T4, mass ratios increase 
with N . In the confined phase, one can compare the ratio 
of the temperature dependent string tension to its value 
at zero temperature. This can be done for N = 2 0, 3 
[To). 4 and 6 [T^. [pi. For example, comparing N = 4 
and 6, Fig. 5 of 13J shows that at fixed T/Td, this ratio 
of string tensions increases as N does. In the deconfincd 
phase, Fig. 7 of the most recent work by Lucini, Teper, 
and Wenger ^| shows that the ratio of the Debye mass 
at T£, to Trf, is significantly larger for N — 8 than for 
N = 3. Thus, the lattice data strongly indicates that the 
N = 00 theory is not exactly at the Gross-Witten point; 
if it were, at finite N one would expect that mass ratios 
would decrease, as T — > Td, in a manner which is nearly 
A^-independcnt . 

A recent, lengthy calculation by Aharony et al. shows 



that on a small sphere, at infinite N the deconfining tran- 
sition is of first order [5]]]. If g is the gauge coupling, 
then as N — > 00, g 2 N is a number of order one, while 
g 2 N <C 1 on a small sphere. In this limit, the coefficient 
of the quartic coupling in the loop potential, l|63|) . is K4 ~ 
[g 2 Nf [3; H3 finds that k 4 « -(l/10)(.g 2 A7(47r)) 2 . 
Although it is difficult to know how to normalize the co- 
efficient, its value suggests that on a small sphere, the 
theory is near the Gross-Witten point. Perhaps this re- 
mains true in an infinite volume. 

With optimism, then, we assume that in the sense of 
Sec. IIII El the deconfining transition at N = 00 is close 
to, but not at, the Gross-Witten point. We suggest that 
at finite N, fluctuations can drive the theory much closer 
to this point: that it is, in the space of all possible cou- 
plings, an infrared stable fixed point. Couplings only flow 
due to fluctuations, however, and in matrix models these 
are manifestly ~ 1/N 2 . Thus this can only happen in a 
region of temperature which shrinks as N increases, for 
\T-T d \/T d ^l/N 2 . 

A comparison of lattice simulations at different values 
of N provides some evidence for a critical window which 
narrows with increasing N. Consider the decrease of the 
string tension near Td- Define a temperature T\m as the 
point at which the string tension is half its value at zero 
temperature, a{T^/ 2 ) = 0.5 c(0). The lattice finds that 
the corresponding reduced temperature, 5t\/2 = {T d — 
T-1/2) /Td, decreases as N increases: 8t\/2 ~ 0.2 at N = 2 
(Fig. 3 of 9]), w 0.08 when N = 3 (Fig. 5 of E3), and 
« 0.05 for N = 4 (Fig. 5 of [H). These three points can 
be fit with 

St 1/2 w 0.8/N 2 . (83) 

This really is a narrow window, with 6ti/ 2 ~ 0.02 for 
N = 6. 

Unfortunately, a loop model cannot predict the ratio of 
the string tension at Td, to zero temperature, because we 
do not know what m 2 corresponds to zero temperature. 
For three colors, this ratio is 0.12 10]. A loop model 
can predict the ratio of the string tension, to the Debye 
mass, at Td- At infinite N , as T — * Td the string tension 
vanishes as a(T) ~ \T - T d \ 1/2 , while the Debye mass 
vanishes more slowly, moebye ~ iT-iy 1 / 4 and O 
0. That a(T)/m Debye (T) ~ \T d - as, T —> T d 

is unlike ordinary second order transitions, where the 
analogous ratio is constant in a mean field approxima- 
tion. For three colors, a{T d /)/m Debve {T+) w 0.93 ± 0.02 
at the N — 'i analogy of the Gross-Witten point, and 
w 0.99 ± 0.03 when n 3 = 0.146, Sec. ITVBl For four 
colors, it is w 0.63 ± 0.02 at the A^ = 4 analogy of the 
Gross-Witten point [50|. and decreases with increasing 
N. 

The most direct way to test if the transition for A^ > 
4 is close to the Gross-Witten point is to compute the 
value of the renormalized fundamental loop (very) near 
Td 0, , and see if it is close to \ . The change in the 
masses near Td is probably more dramatic, though: since 
the latent heat grows as ~ n 2 [HEil, such a decrease, 
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especially in a window which narrows as N increases, 
would be most unexpected. 

To rigorously demonstrate the connection between 
loop models and deconfinement, it will be necessary to 
include fluctuations at finite N. At the beginning of Sec. 
IIIII we blithely ignored the kinetic terms of Sec. Ill CI pro- 
ceeding immediately to a mean field theory, which is an 
expansion in bare quantities. Including fluctuations will 
require the analysis of the renormalized theory. Fluctua- 
tions are most important in the lower critical dimension, 
which is two |35l |49| , and is relevant for the deconfining 
transition in three spacetime dimensions |52j . The effect 
of fluctuations can also be computed by working down 
from the upper critical dimension, which is four. In this 
vein, we comment that the entire discussion of Sec. IIIII is 
dominated by the Vandermonde potential, which arises 
from the measure of the matrix integral. In the contin- 
uum limit, terms in the measure can be eliminated from 
perturbation theory by the appropriate regulator |35j . 



The measure can contribute non-perturbatively, though: 
see, e.g ., Eq. (73) of Billo, Caselle, D'Adda, and Panz- 
eri |3j|. In general, it will be interesting to develop the 
renormalization group for what may be a new universal- 
ity class, about the Gross- Witten point. 
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FIG. 1: The N — oo effective potential at the Gross- Witten 
point. When £n < the potential vanishes identically. In 
Figs. 1, 2, and 3, there is a discontinuity, of third order, at 
In = \, although it is only apparent here. 




FIG. 2: The N = oo effective potential for K4 — — 1, where 
the transition is strongly first order. 
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FIG. 3: The Vandermonde potential, Vvdm, for N = 2 (dot- 
ted line), N = 3 (dashed line) and N — oo (full line). 




FIG. 4: The N = 3 effective potential, at the transition, for 
^3 = (full line) and 0.146 (dashed line). 
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FIG. 5: Expectation values for the 3, 6, 8, and 10 representa- 
tions of SU(3), for K3 = (full lines) and 0.146 (dashed lines), 
as a function of the ratio of m 2 to its critical value. For better 
visibility, the expectation value of the octet has been shifted 
by +0.1 




FIG. 6: Difference loops for 6, 8, and 10 representations of 
517(3), for K3 = (full lines) and 0.146 (dashed lines), as a 
function of the ratio of m 2 to its critical value. 
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FIG. 7: The effective mass squared, for the fundamental loop 
in the deconfined phase, versus its expectation value. The 
three curves are N = oo about the Gross- Witten point, all 
Ki = (full line); N = 3 with all m = (dashed line); and 
N = 3 with k 3 = 0.146 (dotted line). 



